Stable  sets,  corner  polyhedra  and  the  Chvatal  closure 


Manoel  Campelo 

Departamento  de  Estati'stica  e  Matematica  Aplicada, 
Universidade  Federal  do  Ceara,  Brazil, 
mcampelo@lia.ufc.br. 

Gerard  Cornuejols 

Tepper  School  of  Business,  Carnegie  Mellon  University,  Pitsburgh,  PA  15213, 

gcOv@andrew.cmu.edu. 

January  2009,  revised  July  2009. 


Abstract 

We  consider  the  edge  formulation  of  the  stable  set  problem.  We  characterize  its 
corner  polyhedron,  i.e.  the  convex  hull  of  the  points  satisfying  all  the  constraints 
except  the  non-negativity  of  the  basic  variables.  We  show  that  the  non-trivial 
inequalities  necessary  to  describe  this  polyhedron  can  be  derived  from  one  row  of 
the  simplex  tableau  as  fractional  Gornory  cuts.  It  follows  that  the  split  closure  is  not 
stronger  than  the  Chvatal  closure  for  the  edge  relaxation  of  the  stable  set  problem. 
Keywords:  Stable  set,  corner  polyhedron,  Chvatal  closure,  odd  cycle  inequality. 


1  Introduction 

Consider  a  simple  graph  G  =  (V,E),  where  V  and  E  are  the  sets  of  n  vertices  and  m 
edges  of  G,  respectively.  A  stable  set  (independent  set,  vertex  packing)  of  G  is  a  set  of 
pairwise  non-adjacent  vertices.  For  the  sake  of  simplicity,  we  are  going  to  assume  that  G 
has  no  connected  component  defined  by  a  single  vertex.  Then,  a  stable  set  corresponds 
to  an  n- dimensional  binary  vector  x  that  satisfies  xu  +  xv  <  1,  for  all  uv  €  E.  The  set  of 
stable  sets  of  G  can  be  described  by  the  mixed  integer  linear  set 

5(G)G(^)eZ;xR™:  Ax  +  y  =  1},  (1) 

where  A  is  the  edge-vertex  incidence  matrix  of  G,  1  is  a  vector  of  ones,  x  and  y  are 
vectors  of  variables  indexed  by  the  vertices  and  the  edges  of  G,  respectively.  The  convex 
hull  of  S(G)  is  called  the  stable  set  polytope.  The  set  obtained  from  S(G)  by  relaxing 
the  integrality  constraints  will  be  denoted  by  R(G). 

An  important  reason  for  studying  methods  to  handle  contraints  of  type  (1)  is  that  they 
can  model  restrictions  appearing  in  many  optimization  problems.  Indeed,  notice  that  set 
packing  and  partitioning  problems  can  be  transformed  into  vertex-packing  problems  on 
the  intersection  graph. 
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A  general  approach  to  deal  with  mixed  integer  sets  consists  of  solving  the  linear  relax¬ 
ation  and  adding  cutting  planes.  Most  commonly,  these  cutting  planes  are  derived  from 
integrality  arguments  applied  to  a  single  equation.  Intersection  cuts  are  more  general  [4], 
Recently,  cutting  planes  derived  from  two  or  more  rows  have  attracted  renewed  interest 
[2,  6,  9],  and  it  has  been  shown  that  they  can  better  approximate  the  feasible  set  [3,  5]. 
Whether  they  are  derived  from  one  or  more  rows  by  integrality  arguments,  the  known 
cuts  have  in  common  the  feature  that  they  are  valid  for  the  corner  polyhedra  defined  by 
Gomory  [14], 

In  the  remainder,  let  B  stand  for  the  set  of  all  (feasible  or  infeasible)  bases  of  the 
constraint  matrix  [A  /].  The  corner  polyhedron  associated  with  a  basis  B  e  £>,  to  be 
denoted  corner (B),  is  the  convex  hull  of  the  feasible  points  of  the  relaxation  of  (1)  where 
the  non-negativity  constraints  on  the  basic  variables  are  discarded.  A  valid  inequality  for 
corner (B)  is  valid  for  (1). 

In  this  paper,  we  study  the  corner  polyhedra  for  (1)  and,  consequently,  investigate  cuts 
that  can  be  derived  from  a  basic  solution  of  the  linear  relaxation  of  (1).  We  show  that  all 
non-trivial  valid  inequalities  necessary  to  the  description  of  corner (B)  can  be  derived  from 
one  row  of  (1)  as  a  Chvatal- Gomory  cut  [8,  13].  In  addition,  we  relate  the  intersection  of 
the  corner  polyhedra  associated  with  all  bases,  given  by 

corner (B)  =  |^|  corner (R), 
bgb 

to  the  split,  Chvatal  and  {0,  l/2}-Chvatal  closures  relative  to  S(G).  The  Chvatal  closure 
is  the  intersection  of  R(G)  with  all  the  Chvatal  inequalities,  i.e.  inequalities  of  the  form 
[XA\x  <  |_A1J  with  A  <E  [8].  When  A  is  restricted  to  be  in  {0,l/2}m,  the  {0,1/2}- 
Chvatal  inequalities  and  closure  are  similarly  defined  [7]. 

It  follows  from  [1]  that  the  polyhedron  corner(R)  is  included  in  the  split  closure.  In 
turn,  the  split  closure  is  included  in  the  Chvatal  closure,  which  is  itself  contained  in  the 
{0,  l/2}-Chvatal  closure.  Here,  we  show  the  converse  inclusions  and  obtain  the  following 
theorem. 

Theorem  1.  For  the  stable  set  formulation  (1)  ,  the  {0,  l/2}-Chvatal  closure,  the  Chvatal 
closure,  the  split  closure  and  corner(B)  are  all  identical  to 

{(x,y)  >  0  :  Ax  +  y  =  1,  <  - — - — ,V  induced  odd  cycle  C  of  G}. 

vec  2 

In  this  paper,  a  cycle  is  a  sequence  C  =  (r>o,  ei,  iq, . . . ,  e*,,  ty  =  Vo),  where  iq,  tq, ...,  vj~ 
are  distinct  vertices  and  e*  =  for  i  —  1,  2, . . . ,  k,  are  distinct  edges.  Such  a  sequence 

is  often  called  a  circuit  in  the  literature  [11].  The  cycle  C  is  induced  if  every  edge  of  the 
graph  connecting  two  vertices  of  C  is  an  edge  of  G.  It  is  odd  if  the  number  k  of  vertices 
(and  edges)  is  odd.  In  the  remainder,  as  already  used  in  Theorem  1,  |G|  stands  for  the 
number  of  vertices  of  C.  Also,  we  simply  write  v  G  C  or  e  G  C  to  mean  that  a  vertex  v 
or  an  edge  e  is  in  C. 

In  order  to  obtain  Theorem  1,  we  first  describe  a  graph  related  to  the  basic  matrix 
in  the  next  section.  Section  3  is  devoted  to  the  characterization  of  the  corner  polyhedra, 
whereas  Section  4  relates  them  to  the  split,  Chvatal  and  {0,  l/2}-Chvatal  closures. 
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2  The  basic  graph 

Consider  any  basis  B  e  B,  feasible  or  infeasible,  and  the  corresponding  nonbasic  matrix 
N.  Let  xB,yB  and  x^.y^  represent  the  basic  and  nonbasic  variables,  respectively.  The 
vertices  and  edges  indexing  these  variables  will  be  respectively  called  basic  and  nonbasic. 
Let  us  group  the  equations  of  (1)  according  to  the  basic  and  nonbasic  edges  to  get 

Ax  +  yN  =  1,  (2) 

Ax  +  yB  =  1,  (3) 


where  A 


A 

A 


Since  the  non-negativity  of  the  variables  yB  is  discarded  in  the  defini¬ 


tion  of  corner(L>),  constraints  (3)  become  redundant.  So,  let  us  focus  on  the  first  group 
of  constraints.  Toward  this  end,  we  define  the  basic  graph  Gb ,  which  is  obtained  from  G 
by  removing  the  basic  edges.  Notice  that  A  is  the  edge- vertex  incidence  matrix  of  GB. 
Actually,  this  relation  can  still  be  considered  in  the  trivial  case  where  all  edges  are  basic. 

Let  Ci  =  ( Vi,Ei ),  %  =  1,2,...,  A;,  be  the  connected  components  of  GB  that  are  not 
defined  by  a  single  vertex.  If  Ai  is  the  \Et\  x  \Vt\  incidence  matrix  of  Gl,  then  A  can  be 
organized  as 


1 

1 

A2 

1 

1  0 

1 

Ak  ! 

For  every  i  =  1,2 , ,k,  let  us  partition  Ai  =  [Bi  IVj],  where  Bi  and  Nt  comprise  the 
columns  of  Ai  that  are  indexed  by  the  basic  and  nonbasic  vertices,  respectively. 


Lemma  2.  For  every  i  —  1,  2, . . . ,  k,  the  matrix  Bi  is  square  and  invertible. 

,  where  B  and  B  are  submatrices  of  A  and  A,  respec¬ 
tively.  Since  B  is  a  basis,  B  must  be  invertible.  Therefore, 


Proof:  We  have  that  B  = 


B  0 
B  I 


B  = 


Bi 


Bo 


Bu 


In  addition,  the  rows  and  columns  of  each  Bt  must  be  linearly  independent.  So,  each  Bi 
is  square  and  invertible.  □ 


Now,  we  show  that  each  connected  component  of  Gb  is  a  tree  or  a  1-tree.  A  1-tree  is 
a  connected  graph  with  exactly  one  cycle. 

Lemma  3.  For  every  B  €  B,  each  connected  component  of  GB  is  either  a  tree  or  a  1  -tree 
with  an  odd  cycle.  Each  tree  has  exactly  one  nonbasic  vertex.  The  vertices  of  every  1  -tree 
are  all  basic. 
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Proof:  If  a  connected  component  of  GB  is  not  a  tree,  it  must  be  a  component  C*,  for  some 
i  G  {1,  2, . . . ,  k}  such  that  \Ei\  >  \Vi\.  By  Lemma  2,  this  means  that  \Ef  =  \Vf\  and  ZZ*  is 
the  incidence  matrix  of  Ci.  Therefore,  Ct  is  a  tree  together  with  an  extra  edge.  Moreover, 
the  unique  cycle  of  Ci  must  be  odd  because  Bi  is  invertible  [15]. 

The  remaining  part  of  the  statement  follows  by  observing  the  structure  of  B.  It  shows 
that  every  basic  vertex  is  in  some  component  C-i,  which  has  exactly  \V,\  —  \Ei\  nonbasic 
vertices.  □ 

From  now  on,  let  Ib  Q  {1,2, . . . ,  k}  be  the  set  of  indices  of  the  1-tree  components  of 
Gb-  For  i  E  Is,  the  kernel  of  Ci,  denoted  by  «j,  is  its  unique  cycle. 

3  The  Corner  polyhedron 

The  corner  polyhedron  associated  with  a  basis  B  of  the  stable  set  formulation  (1)  is  the 
convex  hull  of  the  points  satisfying  the  relaxation  where  the  non- negativity  constraints 
on  the  basic  variables  are  discarded,  that  is 

corner  (ZZ)  =  conv{(x,  y)  :  Ax  +  y  =  1,  Xn  >  0,  >  0,  x  G  Zn}.  (4) 

Using  (2)— (4) ,  it  follows  that 

corner  (ZZ)  =  conv{(x,  y)  :  Ax  +  =  1,  Xn  >  0,yw  >  0,x  E  TIP,  (5) 

Ax  +  yB  =  1}-  (6) 

Since  ys  is  not  restricted  in  sign  and  appears  only  in  constraints  (6),  we  have  that 

corner (5)  =  {(x,yN,yB)  ■  ( x,yN )  G  convjAx  +  yN  =  1,  xN  >  0,  yN  >  0,  x  G  7Ln}, 

Ax  +  ys  =  1}. 

Notice  that  this  expression  can  also  comprise  the  trivial  case  where  all  edges  are  basic. 
In  this  case,  matrix  A  and  vector  yN  do  not  exist,  which  leads  to  corner (U)  =  {(x,y)  : 
Ax  +  y  —  1,  x  >  0}. 

By  defining  P(B)  =  {x  G  :  Ax  <  1,  Xn  >  0},  we  can  simplify  the  expression  of  the 
corner  polyhedron  as 

corner (B)  =  {(x,y)  :  x  G  conv(P(ZZ)),  Ax  +  y  =  1}.  (7) 

Also,  we  can  rewrite 

P{B)  =  {x£Zn:Ax<  b},  (8) 

where 


'B 

N  ' 

,  b  = 

1 

6 

-21 

0 

and  A  =  [B  N]  is  a  partition  of  A  according  to  (xb,xn). 

Using  expression  (8),  we  can  view  P(B)  as  the  set  of  integral  solutions  of  a  system 
defined  by  the  edge- vertex  incidence  matrix  A  of  a  bidirected  graph  [12].  Such  a  graph  is 
formed  by  GB  together  with  a  loop  for  each  nonbasic  vertex. 
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The  definition  of  cycle  in  graphs  extends  directly  to  bidirected  graphs  [12].  However,  a 
cycle  in  a  bidirected  graph  is  called  odd  if  it  has  an  odd  number  of  edges  whose  endpoints 
are  both  positive  or  both  negative  according  to  the  incidence  matrix  [11,  12]. 

By  Lemma  3,  each  connected  component  of  the  bidirected  graph  associated  with  P(B) 
has  exactly  one  odd  cycle  (a  loop  was  included  in  each  tree  component  of  GB).  This 
property  can  be  used  to  readily  describe  the  convex  hull  of  P(B)  by  means  of  the  following 
result  by  Gerards  and  Schrijver  [12].  See  also  Schrijver  [11,  Section  68.6b],  where  an  odd 
K. 4-sub  division  is  defined  as  a  subdivision  of  the  complete  graph  with  four  vertices  such 
that  each  triangle  has  become  an  odd  cycle. 

Theorem  4  ([12,  11]).  If  M  is  thepx  q  edge-vertex  incidence  matrix  of  a  bidirected  graph 
with  no  odd  AV subdivision  and  b  G  IP ,  then 

conv{x  G  7Lq  :  Mx  <  b}  = 

{x  G  M9  :  Mx  <  b,  ^  EE  Mev%v  —  [_2  E  be J,  for  each  odd  cycle  C}. 

v&C  esc  esc 

Since  the  bidirected  graph  defining  P(B)  is  such  that  each  connected  component  has 
exactly  one  cycle,  it  does  not  contain  a  A^-sub division.  Therefore,  the  above  theorem 
gives  the  following  statement. 

Corollary  5.  For  every  B  G  B,  conv(P(B ))  =  (x  :  Ax  <  1,  Xn  >  0,  < 

(|k*|  -  l)/2,  Vi  G  Is}. 

Corollary  5  together  with  (7)  imply  the  following  characterization  of  the  corner  poly¬ 
hedron. 

Theorem  6.  For  every  B  G  B,  the  corner  polyhedron  of  (1)  associated  with  B  is 
corner(B)  =  <  (x,  y)  :  Ax  +  y  —  1,  xjv  >  0,  vn  >  0,  E  xv  <  <  1,  Vi  G  IB 

I  V &Ki 

Theorem  6  can  also  be  proved  without  appealing  to  Theorem  4  by  characterizing 
the  matrix  B~1,  which  essentially  gives  the  basis  inverse.  Actually,  we  can  show  that 
each  column  of  the  coefficient  matrices  of  B~1Ax  +  B~1y^  =  B~ll  has  its  entries  either 
in  {0,±1}  or  in  {0,  ±1/2}.  This  special  structure  also  leads  to  the  characterization  of 
conv(P(B))  and,  consequently,  the  statement  of  Theorem  6. 

Remark  7.  Each  odd  cycle  inequality  appearing  in  the  description  of  the  corner  polyhe¬ 
dron  is  a  (0, 1/2 }-Chvatal  inequality.  Indeed,  it  can  be  obtained  with  the  Chvatal  procedure 
by  taking  a  linear  combination  of  the  constraints  xu  +  xv  <  1,  for  all  e  =  uv  G  Ki,  with 
multipliers  \e  =  1/2,  for  all  e  G  Kt. 

The  above  results  imply  that  each  non-trivial  inequality  necessary  to  describe  the 
corner  polyhedron  is  a  fractional  Gomory  cut  that  can  be  derived  from  a  row  of  the 
simplex  tableau  related  to  a  basic  vertex  in  the  cycle  of  a  1-tree  component.  Even  if  we 
consider  the  tighter  relaxation  of  (1)  where  only  constraints  yB  >  0  are  discarded,  we 
cannot  get  stronger  cuts.  Indeed,  if  we  keep  the  non-negativity  of  the  (basic)  variables 
associated  with  the  vertices,  we  can  still  use  Theorem  4  to  see  that  only  the  trivial 
inequalities  xB  >  0  can  be  obtained  in  addition  to  those  describing  corner (B). 
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4  The  Chvatal  closure 


Theorem  6  and  Remark  7  show  that  the  only  non-trivial  inequalities  needed  to  define 
corner(R)  are  {0,  l/2}-Chvatal  inequalities.  They  all  have  the  form  Y2vec  x,  <  (|C'|-l)/2, 
where  C  is  an  odd  cycle  of  G.  If  this  cycle  is  not  induced  in  G,  the  corresponding  inequality 
is  dominated  by 


vec1 


xv  <  (|C"|  —  l)/2,  where  C'  is  an  odd  cycle  of  G  induced  by  a 
subset  of  vertices  of  C .  This  is  the  key  point  to  characterize  the  Chvatal  closure  of  S(G), 
as  follows. 

Let  C  denote  the  set  of  all  the  induced  odd  cycles  of  G.  For  every  G  E  C,  form  the 
submatrix  Be  of  [A  I]  given  by 


Bc  = 


Ac  0 

0  / 


where  Ae  is  the  edge-vertex  incidence  matrix  of  C  and  /  is  related  to  the  edges  of  G  not 
in  C.  Notice  that  Be  is  invertible  and  Bffl  e  {1/2,  l}m,  which  implies  that  Be  is  a 
feasible  basis.  Let  us  denote  by  B+  =  {B  E  B  :  L>-11  >  0}  the  set  of  feasible  basis  and 
by  Be  =  {Bc  :  C  E  C}  its  subset  corresponding  to  C.  Also,  for  any  B'  C  B,  let  corner^') 
stand  for  the  intersection  of  the  corner  polyhedra  associated  to  all  bases  in  £>',  that  is, 


corner (15/)  =  f^|  corner(R). 


Lemma  8.  corner{B)  C  corner(B+)  C 


B&B' 


x 


Proof:  We  only  need  to  prove  the  second  inclusion.  First,  note  that  corner (B+)  C 
corner (/)  C  x  Mm.  Now,  for  each  edge  e  =  uv  of  G,  consider  the  submatrix  Be 
of  [A  I]  defined  by  the  column  of  A  indexed  by  u  and  the  columns  of  /  indexed  by 
E  \  {e}.  Notice  that  Be  E  B+.  In  addition,  by  Theorem  6,  ye  >  0  is  a  valid  inequality  for 
corner (Re).  Since  corner (R+)  C  nee^corner(Zje),  the  result  follows.  □ 


\C\ 


Lemma  9.  corner(B)  =  corner(B+)  =  S(G),  where 

S(G)  =  <  (x,  y)  >  0  :  Ax  +  y  =  1,  ^  xv  < 

[  vec 

Proof:  If  C  —  0,  the  result  trivially  follows  by  Theorem  6  and  Lemma  8.  The  same 
statements  imply  that 


-yc  e  c 


corner  (7?)  C  corner  (£>+)  C  corner  (R^)  fl 


)  =  S(G). 


To  obtain  equalities  in  the  expression  above,  it  suffices  to  note  that  there  is  always  a  subset 
of  the  vertices  of  any  non-induced  odd  cycle  C  that  defines  an  induced  odd  cycle  C'  E  C. 
Then,  the  inequality  X^ueC'  Xv  —  (1^1  —  1) /2  is  tighter  than  'f2vGCxv  A  (|C|  — 1)/2,  which 
shows  that  S(G)  C  corner (B).  □ 

Lemma  9  and  the  fact  that  all  inequalities  needed  in  the  description  of  corner (£>)  are 
{0, 1/2 {-Chvatal  inequalities  lead  to  the  following  more  complete  version  of  the  theorem 
stated  in  the  introduction. 
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Theorem  10.  For  the  stable  set  formulation  (1),  S(G),  the  {0, 1/2 }-Chvatal  closure,  the 
Chvatal  closure,  the  split  closure,  corner{B+)  and  corner(B)  are  identical. 

We  should  stress  that  the  equivalences  stated  in  Theorem  10  may  not  be  valid  for 
other  formulations  of  the  stable  set  problem.  For  instance,  for  a  quasi-line  graph  the  split 
closure  of  the  clique  formulation  of  the  stable  set  problem  coincides  with  the  convex  hull 
of  the  integer  solutions  and  is  strictly  contained  in  the  Chvatal  closure  [10]. 
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